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using  target  values  as  model  inputs  for  not  only  the  beqinning  of  a  conflict, 
but  also  for  times  extending  into  the  conduct  of  a  conflict.  New  technoloqy 
has  been  applied  to  the  problem  and  the  success  achieved  is  reviewed. 
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I.  THE  PROBLEM 

Military  Planning  for  the  United  States  primarily  involves  preparation 
for  an  aggressive  response  to  an  attack  made  against  the  United  States  or 
allied  forces.  This  basically  defensive  posture  all  but  denies  the 
United  States  the  advantages  inherent  in  surprise,  and  requires,  as  a 
counter  to  the  initiative  granted  an  attacker,  that  the  United  States 
maintain  operational  effectiveness  and  vigilance  second  to  none. 

To  achieve  maximum  military  effectiveness  given  fiscal  constraints 
requires  careful  consideration  in  planning  and  budgeting.  Men  without 
weapons  will  neither  intimidate  nor  stop  an  enemy,  neither  will  men  with 
weapons  but  without  munitions  or  a  means  of  maneuver  (transport).  As  an 
example,  the  United  States  Air  Force  must  weigh  each  dollar  spent  to  see 
if  it  should  be  best  utilized  for  personnel  (acquisition,  retention,  and 
training),  aircraft  (acquisition  and  maintenance),  fuel,  or  munitions  in 
order  to  best  fulfill  its  many  missions. 

The  research  described  herein  addresses  the  Air  Force  nonnuclear 
munitions  procurement  model.  This  model  requires  as  input: 

1.  A  target  list  to  include  target  military  values,  defensive  capabil- 
ities, and  factors  describing  the  potential  for  confirming  tarqet 
kill; 

2.  Numbers  of  friendly  aircraft  and  munitions  available; 

3.  Probabilities  of  encountering  differing  weather  conditions;  and 

4.  The  effectiveness  of  specific  munitions  against  various  targets 
when  delivered  in  a  weather  condition  by  a  specified  aircraft. 

Given  these  inputs,  the  model  seeks  to  maximize  the  damage  done  to 

the  enemy  by  planning  the  use  of  the  most  effective  munitions. 

3 


The  model  is  played  separately  for  each  major  Air  Force  theatre. 
Individual  theatre  results  comprise  the  single  largest  input  for  non- 
nuclear  weapons  procurement  decision-making  for  the  Air  Force.  The 
amount  of  money  spent  with  the  aid  of  this  model  is  currently  in  excess 
of  one  billion  dollars  a  year.  Not  only  is  the  model  central  to  Air 
Force  budgeting  and  planning,  but  it  is  also  relied  upon  by  the  major  theatre 
commanders  for  insight  into  scenarios  involving  their  current  missions. 

With  so  much  at  stake  in  its  use,  an  examination  is  in  order  to 
better  understand  the  assumptions,  formulation,  processing,  mathematical 
solution,  and  solution  report  generation  for  this  model  [Ref.  1], 


II.   DESCRIPTION  OF  THE  HEAVY  ATTACK  (HA)  MODEL 

Heavy  Attack  (HA)  is  a  program  which  uses  an  imbedded  nonlinear 
optimizer  to  identify  a  set  of  sortie  allocations  which  maximizes  the 
military  worth  of  targets  killed.   It  has  been  used  since  1974  as  part  of 
a  set  of  computer  programs  that  have  been  known  as  Saber  Mix  or,  more 
recently,  as  the  Nonnuclear  Armament  Plan  (NAP)  models.  These  models 
together  attempt  to  provide  an  optimal  munitions  mix  for  a  given 
specification  of  available  sorties,  targets,  and  other  factors. 

The  HA  model  can  be  viewed  as  consisting  of  an  internal  (optimization) 
model  and  an  external  model  consisting  of  input  and  the  output  sequences. 
This  view  will  be  adopted  in  describing  the  model  and  its  variations 
primarily  because  the  optimization  is  still  modeled,  if  not  processed, 
exactly  as  it  first  was  eight  years  ago,  while  input  and  output  models 
have  changed  (grown)  continuously.  Careful  attention  will  be  paid  to  the 
optimization  model  while  the  input  and  output  models  will  be  examined 
only  to  achieve  an  understanding  of  the  consequences  of  model  input 
aggregation  and  output  unraveling  for  solution  interpretation. 

A.   DESCRIPTION  OF  THE  HA  OPTIMIZATION  FORMULATION  (INTERNAL  MODEL) 

HA  solves  a  sequence  of  internal  models.  Each  of  these  models  has 
the  same  mathematical  structure  (formulation),  and  each  is  a  nonlinear 
optimization  problem  by  virtue  of  its  objective  function.  The  internal 
optimization  model  is  formulated  as  follows: 
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maximize: 


subject  to: 


Z     vj(lW  -  Dj) 


J 

< 


£      Sid   =  Sis    i   =  1,   2,    ...,    I; 


£j  £  Kj  £  Tj>     j   =  1,   2,    ...,   J; 


6 
m 


(Z    (i-  V  s,  d  -  Z    (-V  Si  -)<o, 


m  =   1,    2,    . . .   M; 


Sij  >  °- 


The  objective  function,  which   is  convex  with  respect  to  S.  . 
[Ref.   2:     pp.   8,   9],   quantifies  the  value  of  all   targets  killed   in   a 
time  period  as  follows: 

Z  VVV  -  Dj>- 

where, 

i  =  aircraft  type  index  (i  =  1,  ...,  I); 

j  =  target  type  index  (j  =  1,  ...,  J); 

S • •  =  number  of  sorties,  the  independent  variables; 

V.  =  value  (military  worth)  of  target  type  j; 

0.  =  cumulative  number  of  targets  killed  in  prior  time  periods;  and 

J 

K.^S  —  }  =  number  of  kills  of  target  type  j,  a  nonlinear  function  of  S-- . 

K.  is  defined  to  be  a  function  of  the  independent  variables  S-- : 
j  j 

T.  C.        I 

J.  M   -   avr,  I 1 


Kj<S1j>=^(l-exp(  -T^VS/idV'' 


where, 


T.  =  number  of  type  j  targets; 

J 
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C  =  target  kill  conf  irmabil  ity  parameter  (0  <_  C-  <_  1)  controlling  the 
J    extent  to  which  the  law  of  diminishing  productivity  (as  described 
later)  applies; 

T.         C. 

0.  -  t^-   Tog  (1  -  yj-  0-:);  a  term  (to  be  explained)  added  for  mathematical 
J     j         j     convenience;  and 

P--  =  expected  number  of  type  j  targets  killed  per  type  i  aircraft 
sortie  when  no  other  targets  of  type  j  have  been  previously 
killed  and  when  conditions  of  kill  conf irmabil ity  are  perfect. 

HA  problem  constraints  are  of  three  forms: 

1.  Sorties  available: 

J 
g  S..   =  Sif  i  =  1,  2,  ...,  I; 

where  S.  =  the  number  of  sorties  for  aircraft  type  i. 

2.  Target: 

*j  1  Kj  <  Tj,  j  ■  1,  2,  ....  J; 

where  i.   =  the  lower  bound  on  targets  of  tyoe  j  which  must  be 
killed. J 

3.  Fl ight  composition  with  the  general  form: 


5m 


(Z      (1-  em)  Si   +  L  (-em)  Si   )<0,  b-1.2,  ....  N; 


m 
where, 

m  =  flight  composition  constraint  index  (m  =  1,  2,  ...,  M) ; 

i  =  aircraft  type; 

<s  =  +1,  maximum  (-1,  minimum); 

Jm  =  set  of  targets  for  which  a  maximum  (minimum)  flight 
composition  is  required;  and 

9  =  maximum  (minimum)  proportion  of  sorties  flown  by  aircraft 
type  i  against  targets  included  in  set  J 
[Ref.  2:  pp.  5,  8,  10].  m 

K.  is  computed  assuming  "a  law  of  diminishing  marginal  Droductivity  .  . 

i.e.,  the  number  of  targets  per  sortie  decreases  for  each  successive  sortie' 
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[Ref.  2:  p.  1].  This  underlying  assumption  is  based  on  a  belief  that 
targets  will  be  harder  to  find,  harder  to  surprise,  and  harder  to  kill  as 
the  battle  continues. 

By  restating  target  constraint  inequalities  to  solve  for  z  P^-S--, 
the  problem  can  be  formulated  (see  Appendix  A)  with  linear  constraints 
[Ref.  2:  pp.  6-7,  10-12]  of  the  form: 


i=i   IJ  1J 


where, 


-T. 
Rj  =-Ti  dog  (1  -  Cj))  -  aj;  and, 

-T.  C. 

Rj  =-7^  (1o9  (1  "  *i  J1   ))  "  ai- 

Further  examination  of  K.  provides  significant  insight  into  the 
nature  of  the  function  which  is  being  maximized.   If 

T.  C.       I 


Kj'=^(1"exp("T7(aj  +  SPi^Si^)))' 

then  noting, 


C.  C. 

exp  (-  J-  o.)  »  1  -  -J-  D., 

j  J        j  J 

and  (for  convenience  of  exposition)  letting, 

I 
X.  =  2  P.  .S.  ., 

the  term  K.  -  D.  (from  the  objective  function)  can  be  restated  (see 

Appendix  B)  as  follows: 

T.  C. 

Kj  -  Dj  =  (  ^-  Dj)(l  -  exp  (Xj  (-^L  ))}. 
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The  function  K -  -  D -  is  a  composition  of  two  embedded  functions. 

The  constant  parameter  C-  can  be  used  to  select  either  component.  When 

Cj  is  equal  to  one, 

-X 


'j-DjMTj-DjHl-exp-a). 


When  C.  is  zero, 

Regardless  of  C,  if  X-/T-  is  small  (a  target-rich  environment  exists), 

vJ  \J     vJ 

then  K.  -  D.  =  X,. 

J        J        J 

C.  can  be  interpreted  to  be  a  coefficient  of  conf irmabil ity.  That 
is,  if  conditions  allow  a  pilot  to  confirm  the  effect  of  his  first 
ordnance  drop  before  delivering  any  others  then  C.  should  be  zero.  This 
case  can  be  called  shoot-look-shoot. 

If  conditions  are  such  that  a  pilot  must  deliver  his  ordnance  without 
reference  to  the  success  or  failure  of  any  weapons  that  might  have  been 
delivered  earlier,  then  C.   should  be  one.  This  case  could  be  called 

J 

dump-al 1-ordnance. 

While  the  composition  of  K.  -  D-  can  be  understood  for  the  boundary 
conditions  where  C.  is  equal  to  one  or  zero,  it  is  not  at  all  certain  how 

vJ 

K.  -  D-  should  be  interpreted  in  cases  where  C-  lies  between  the  two 
boundaries,  except  that  it  will  behave  as  some  mixture  of  the  two  embedded 
functions. 

Summarizing,  HA  solves  a  problem  with  a  nonlinear,  but  convex, 
objective  function  subject  to  three  sets  of  linear  constraints.  Function- 
ally it  selects  sorties,  not  weapons 3   to  inflict  the  greatest  damage  upon 
an  enemy.  Constraints  demand  each  aircraft  be  utilized  when  available, 
that  a  specified  range  of  each  target  type  must  be  killed,  and  that 
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specified  sets  of  targets,  associated  with  an  aircraft,  must  be  attacked  by  a 
minimum  (or  maximum)  percentage  of  the  sorties  available  for  that  aircraft. 

B.   DESCRIPTION  OF  THE  JUNE  1982  HA  EXTERNAL  MODEL 

Because  the  optimization  model  is  formulated  in  terms  of  sorties 

(aircraft/target  combinations)  while  the  purpose  of  HA  is  to  provide 

preferred  weapon   information,  the  HA  external  model  largely  concerns 

itself  with  building  sortie  information  (P-.'s)  from  inputs  and  extracting 

from  sortie  solutions  (S-.'s)  information  regarding  preferred  weapons. 

The  HA  model,  external  to  the  optimization  model,  has  been  modified 

repeatedly  in  the  years  since  1974  for  a  variety  of  reasons.  Determining 

the  history  of  these  changes  might  prove  interesting,  but  what  is  more 

important  is  to  understand  HA  in  its  present  form.   (Inputs,  processes, 

and  outputs  for  the  original  HA  model,  as  described  in  the  seminal  paper 

for  the  model  [_Ref.  2:  p.  2],  are  described  in  Appendix  C.) 

Preceding  HA  in  the  execution  of  the  Nonnuclear  Armament  Plan  (NAP) 

model  programs  are  Weaponeer,  Survivor,  and  Selector  [Refs.  3:  p.  5  and 

Ref.  4:  p.  2]. 

"The  first  model,  Weaponeer,  computes  the  expected  target  kills 
per  pass  for  various  aircraft/munition/target  combinations.  .  .  . 
The  second  model,  Survivor,  computes  the  attrition  of  delivery 
aircraft  on  an  iterative  basis  and  determines  the  expected  kills 
per  sortie  ....  The  third  model,  Selector,  selects  the  preferred 
long-list  weapon  for  each  aircraft/target/delivery  band  [weather 
type]  combination.  ..."  [Ref.  4:  p.  2] 

The  models  preceding  HA  in  the  NAP  sequence  are  of  current  interest 

only  insofar  as  they  provide  inputs  to  the  HA  optimization.  What  should 

be  noted,  however,  is  the  concern  for  the  passage  of  time  implicit  in  the 

short  descriptions  of  Weaponeer  and  Survivor  given  above  (and  in  the 

prior  discussion  of  diminishing  marginal  productivity),  as  well  as  the  use 
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of  weather  types  in  Selector.  While  Clasen,  Graves,  and  Lu  [Ref.  2] 
never  discuss  use  of  the  model  over  multiple  time  periods  or  delivery 
(weather)  types,  these  two  elements  are  currently  incorporated  in  the  HA 
model.  (Multiple  time  period  use  is  implied  in  [Ref.  2]  by  inclusion  of 
the  term  0.  in  the  objective  function  and  by  the  adjective  "interval" 
describing  model  inputs  and  outputs,  however,  no  "iteration-over-time" 
scheme  is  described.  The  effects  of  weather  on  sortie  effectiveness  are 
never  mentioned  in  [Ref.  2].) 

Weather  conditions  have  a  profound  impact  on  aircraft/weapon  effective- 
ness. The  HA  external  model  expresses  weather  condition  as  six  discrete 
"bands,"  associated  with  extremely  poor  to  essentially  unlimited  flying 
visibility.  The  June  1982  HA  input  model  builds  sorties  and  associated 
expected  kills  using  a  weighted  average  based  on  expectation  of  weather 
types  in  the  theatre  being  modeled.  The  P..  associated  with  each  S--  is 

J  J 

formed  as  follows: 

numw 


p..  -  r  p  r  p.. i 

"•J   Ti   w  k   ljkw 


where, 

k     =  set  of  all  weapon  types; 

Pw    =  probability  of  a  particular  weather  type; 

^iikw  =  exPected  number  of  kills  by  aircraft  i  loaded  with  weapon  k 
against  target  j  in  weather  condition  w;  and 

numw  =  number  of  weather  condition  types  modeled. 
The  effect  of  this  averaging  prior  to  the  optimization  model  is  two- 
fold. Sorties  are  not  limited  to  a  single  weapon  type;  that  is,  sorties 
are  evaluated  as  aircraft,  not  as  aircraft/weapon  combinations.  More 
importantly,  an  average  quantity  representing  the  aircraft's  effectiveness 
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is  being  sent  to  the  optimizer.  The  result  is  a  selection  of  extremal 
averages  rather  than  the  averaging  of  selected  extremal  values. 

This  can  lead  to  paradoxical  recommendations.  Consider  a  paradigm 
with  two  equally  probable  weather  conditions,  two  aircraft  types,  and  a 
single  target  type.  Given  one  aircraft  is  moderately  effective  against 
the  target  in  both  weather  types  while  the  second  aircraft  is  near 
certain  to  kill  the  target  in  one  weather  type  and  as  certain  to  miss  it 
in  the  other,  the  first  aircraft  is  likely  to  be  selected  by  the  model 
(to  be  flown  in  both  weather  types)  despite  the  dominance  of  the  other 
aircraft  in  the  second  weather  type. 

In  HA  the  effects  of  weighted  averaging  may  be  subtle.  A  situation 
is  likely  in  which  for  four  of  six  weather  types,  the  effectiveness  of  an 
aircraft/best  weapon  combination  against  a  target  type  is  quite  high 
while  in  the  other  two  weather  types,  the  effectiveness  of  this  aircraft 
is  negligible.   If  another  aircraft  has  mediocre  success  in  all  weather 
types  against  the  same  target,  it  might  well  be  selected  by  the  optimizer 
because  the  optimizer  is  not  allowed  to  see  the  "extreme"  effectiveness 
of  the  first  aircraft  against  the  target  type  in  particular  weather 
conditions.  A  weighted  average  model  will  not  necessarily  provide  valid 
answers  for  any  HA  scenario. 

A  second  major  function  of  the  external  model  is  to  accommodate 
processing  over  multiple  time  periods.  The  number  of  time  periods  is 
limited  to  seven  in  the  June  1982  model,  and  each  time  period  is  typically 
defined  to  last  fourteen  to  thirty  days.  Model  parameters  do  change  over 
time.  The  number  of  dead  targets  (D,)  is  accumulated  following  optimization 
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for  each  time  period.  Other  parameters  (in  particular,  target  values) 
are  reinitialized  at  specified  time  intervals. 

The  use  of  multiple  time  periods  in  HA  is  myopic.  That  is,  optimiza- 
tion proceeds  forward  in  time,  in  one  pass,  with  no  backtracking.  (Each 
optimization  seeks  to  inflict  maximum  damage  upon  enemy  targets  in  that 
time  period   without  regard  for  the  effect  "decisions"  made  in  the  current 
time  period  may  have  on  outcomes  for  future  time  periods.) 

It  is  apparent  that  a  target  type  with  greater  than  average  defenses 
will  generally  have  a  higher  than  average  value.  Having  a  higher  target 
value  will  increase  the  likelihood  the  target  will  be  attacked — at 
least  in  a  short  duration  "conflict"--since  we  are  simply  trying  to 
inflict  the  greatest  damage  upon  the  enemy.  However,  by  avoiding  this 
target  type  in  a  longer  duration  conflict  and  killing  other  less  valuable 
and  less  heavily  defended  targets,  the  attacking  aircraft  might  be  used 
to  inflict  more  damage  on  the  enemy,  over  time,  by  surviving  longer  and 
killing  a  large  number  of  less  valuable  targets. 

Lack  of  logistical  constraints   poses  the  most  significant  difficulty 
of  a  practical  nature  for  the  June  1982  model.  HA  uses  a  preferred 
weapon  without  regard  to  its  actual  availability.  The  internal  model 
never  "sees"  weapons--it  only  "sees"  sorties.  During  a  time  interval, 
the  model  will  continue  to  "fly"  a  sortie  type  as  long  as  the  objective 
function  value  is  improved  (subject  to  the  model  constraints).  Even  when 
it  is  evident  the  supply  of  a  preferred  weapon  has  been  exhausted  in 
previous  time  periods,  there  is  no  way  to  prevent  the  June  1982  model 
from  continuing  to  use  that  preferred  weapon  in  subsequent  solutions. 
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C.   DESCRIPTION  OF  THE  JUNE  1982  HA  PROGRAM  CODE 

Appendix  D  is  a  hierarchy  chart  of  the  June  1982  HA  programs.  Sub- 
program NONLIN  and  all  those  programs  strictly  subordinate  to  it  on  that 
chart  comprise  the  nonlinear  optimizer,  while  the  other  programs  manage 
input  and/or  output. 

Calling  disciplines,  variable  naming  conventions,  and  documentation  of 
the  HA  code  are  neither  standardized  nor  consistent.  Variable  names  in 
the  system  frequently  do  not  agree  with  current  system  documentation  or 
with  [Ref.  2],  nor  are  they  consistent  among  subroutines.  (Variable 
names  within  the  optimization  code  are,  for  the  most  part,  consistent 
with  [Ref.  2].)  Appendix  E  details  the  name  changes  of  some  of  the 
variables  pertinent  to  the  interface  between  the  main  program  and  the 
optimizer.  Appendix  F  [Ref.  4:  pp.  4-5]  provides  a  data  flow  for  the 
entire  NAP  process  as  it  was  executed  in  June  1982. 
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III.   DISCUSSION  OF  TARGET  VALUES  IN  HEAVY  ATTACK 

A.   CRITICISM  OF  THE  USE  OF  TARGET  VALUES 

In  a  recent  analysis  of  HA  by  G.  Jenkins  [Ref.  5]  (an  Air  Force 

civilian  employee  and  user  of  HA)  the  manner  in  which  inputs  are  prepared 

for  HA,  as  well  as  the  optimization  technique  used  in  the  model,  are 

discussed.  He  concludes  the  optimization  is  straightforward  and  infers 

all  inputs  except  one  are  derived  objectively  and  correctly.  The  one 

input  he  expresses  concern  over  is  target  values. 

"...  the  entire  process  is  based  on  optimization  of  firepower 
scores.  This  model  calls  it  military  worth  which  is  probably  more 
correct,  since  command  centers  and  runways  really  don't  fit  in  the 
context  of  firepower  scores  as  do  tanks.  Nevertheless  the  purpose 
of  this  entire  process,  which  is  pursued  with  meticulous  objectivity 
throughout  each  set  .  .  .  boils  down  to  optimizing  the  relative 
subjective  worth  of  target  values.  Granted,  this  methodology  is 
probably  more  credible  than  the  proverbial  smoke-f il led-room 
approach;  however,  it  indicates  that  there  is  still  room  for 
improvement"  [Ref.  5:  p.  12] 

He  continues, 

"...  there  is  room  for  concern  over  whether  the  model's  solution 
is  credible  at  all.  Maybe  the  model  merely  serves  as  a  guide  to 
quantify  and  substantiate  some  decision-maker's  intuitive  feeling, 
so  that  he  may  proceed  with  what  he  always  wanted  or  'knew'  to  be 
true.  Or,  maybe  there  is  a  sincere  interest  in  gaining  insights 
into  the  combat  process.   If  so,  there  is  plenty  of  room  for 
improvement  in  this  model.  .  .  .  "[Ref.  5:   p.  13] 

Jenkins  is  not  alone  in  expressing  concern  over  the  use  of  target 
values  in  computing  an  optimal  munitions  mix  for  air  forces. 

NATO  first  implemented  the  entire  Sabre  Mix  (now  called  NAP)  method- 
ology with  few  alterations,  but  later  removed  target  values  from  the 
formulation.  The  new  formulation  simply  seeks  to  maximize  the  total 
number  ot  targets  killed  subject  to  a  (new)  constraint  for  each  target 
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type  mandating  predetermined  target  type  kill  proportions  relative  to 

other  target  types  (a  proportionality  constraint).  The  new  objective 

function,  new  constraints,  and  other  problem  constraints  (much  as 

originally  found  in  HA)  are  optimized  with  a  linear  programming  code 

[Ref.  6:  pp.  4-11].  The  reasons  for  embarking  on  the  reformulation  were 

given  as  follows: 

".  .  .  While  it  is  recognized  that  the  concept  of  military  worth 
plays  a  central  role  in  tactical  air  mission  planning  models  with 
short  planning  periods,  this  approach  was  abandoned  for  the  following 
reasons: 

1.  Considerable  difficulty  arose  from  trying  to  assign  credible 
military  worth  functions  to  the  different  target  types. 

2.  It  was  considered  advantageous  to  let  proportions  in  which 
targets  of  different  target  types  are  killed  constitute  input 
to,  and  not  output  from,  a  model  designed  to  contribute  to 
the  solution  of  a  logistical  problem."  [Ref.  6:  p.  3] 

The  concerns  expressed  regarding  the  use  of  target  values  in  HA  can 
be  summarized  as  follows.  Target  values  are  subjectively  derived  and, 
perhaps,  invite  manipulation  of  the  model  by  those  executing  it.  They 
are  hard  to  assign  in  a  credible  fashion.  They  remain  the  only  subjective 
input  in  an  otherwise  objective  process,  and  their  use  might  confuse  the 
distinction  between  tactical  and  logistical  decision  making. 

Discussions  with  Major  F.  Cooper,  the  officer  currently  charged 
with  executing  the  NAP  models  [Ref.  7],  reveals  another  perspective  on 
the  use  of  target  values  in  HA.  The  HA  model  is  run  separately  for  each 
major  theatre  (Europe,  Pacific,  S.E.  Asia,  etc.)  by  a  theatre  project 
team  consisting  of  modeling  personnel  (from  Major  Cooper's  office, 
AF/XOX)  and  military  contingency  planners  and  intelligence  experts 
currently  assigned  to  the  command  responsible  for  the  theatre  being 
examined.  Projected  targets  are  identified  and  target  values  are 
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assigned  for  the  first  time  in  a  given  scenario.  Targets  are  evaluated 
using  a  scale  of  zero  to  twenty  (relative  to  a  tank  platoon  which  is 
given  a  base  score  of  one).  The  model  is  then  executed  for  one  time 
period.  Based  on  the  results  for  that  first  time  period  and  an  examina- 
tion of  the  projected  replenishment  capabilities  of  the  allied  and 
opposing  forces,  target  values  are  assigned  for  the  second  time  period. 
This  process  is  repeated  until  the  model  has  been  run  for  all  time 
periods.  This  cycle — process,  evaluate,  modify,  and  reprocess — has 
historically  developed  credibility  for  the  model  among  its  theatre  users. 
Manual  intervention  also  allows  target  values  to  be  manipulated  so  as  to 
discourage  inadmissable  solutions  produced  by  the  June  1982  HA  model 
(because  of  missing  constraints).  Such  manual,  judgmental  manipulation 
is  done  with  the  full  knowledge  and  concurrence  of  the  theatre  project 
team,  and  its  impact  on  model  realism  is  carefully  evaluated  before 
continuing  with  the  next  time  period.  This  careful  evaluation  of  results 
for  each  time  period  often  requires  twenty  or  more  model  runs  for  a  given 
scenario.  However,  the  model  is  thus  not  allowed  to  use  weapons  no  lonqer 
procurable,  or  to  use  weapons  in  quantities  greater  than  can  be  procured, 
or  to  fly  an  aircraft  to  targets  outside  the  aircraft's  range.  When  a 
theatre  study  is  completed,  the  results  reflect  the  project  management 
team's  consensus  and  its  total  combined  military  judgment. 

B.   A  CASE  IN  BEHALF  OF  THE  USE  OF  TARGET  VALUES  IN  HA 

The  use  of  target  values  in  HA  is  viewed  in  vastly  different  ways 
by  knowledgeable  people  familiar  with  the  model.   In  reviewing  the 
methodologies  used  by  USAF  nonnuclear  munitions  mix  models,  Jenkins  found 
the  use  of  target  values  in  HA  to  be  a  weak  link  in  an  otherwise  strong 

22 


chain.  Loritzen  finds  their  use  difficult  and  inappropriate,  while 

Cooper  finds  them  absolutely  necessary  if  the  model  is  to  reflect  both 

the  commander's  priorities  (as  understood  by  his  contingency  planners  and 

intelligence  experts)  and  real-world  constraints  not  accommodated  by  the 

model's  formulation  [Refs.  5,  6,  7], 

Jenkins  infers  HA  should  replace  the  use  of  target  values  with  a 

two-sided  game  reflecting  optimal  strategies, 

".  .  .  we  can  see  that  this  model  is  not  totally  unlike  the  other 
models  ....  However,  there  is  no  direct  gaminq  structure.  This 
model  is  purely  a  one-sided  affair  and  the  only  service  provided  by 
the  opponent  is  the  supply  of  targets  and  an  unaffected  attrition 
rate.  The  staging  of  the  scenario  and  reconstitution  features  of 
the  model  could  involve  certain  optimal  strategy  games,  but  there 
is  no  evidence  of  this  in  the  process  .  .  ."  [Ref.  5:  p.  12] 

While  the  above  is  intriguing,  developing  the  two-sided  game  preferred 
by  Jenkins  would  require  a  scale,  a  set  of  values  with  which  both  sides 
will  be  originally  endowed  and  by  which  each  side  will  either  gain  or 
lose,  depending  on  their  strategies  and  initial  endowments.  It  is 
apparent  that  this  scale  of  values  will  have  to  relate  to  target  values. 
If  it  is  granted  that  values  for  potential  U.S.  targets  are  not  derived 
easily,  then  how  much  harder  must  it  be  to  evaluate  targets  on  a  scale 
applying  to  both  the  U.S.  and  its  opponent(s)? 

NATO  removed  target  values  from  the  formulation  for  its  strategic 
aircraft  munitions  mix  model,  but  then  apparently  had  a  difficult  time 
validating  the  proportionality  constraint  used  in  the  new  formulation. 
Studying  the  linear  programming  dual  of  their  formulation  reveals  the 
dual  variables  are  "marginal  implied  military  worths  of  the  targets" 
[Ref.  6:  p.  15J.  In  fact,  Loritzen  concludes  that  in  order  to  get  valid 
results  from  the  new  formulation  "the  adjustment  of  proportionality 
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factors  will  necessarily  have  to  be  done  manually  by  military  experts 
until  a  situation  is  reached  where  the  implied  military  worths  do  not 
differ  significantly  from  their  estimates  ..."  [Ref.  6:  p.  15]. 

If  reformulating  HA  as  a  game  or  reformulating  it  using  target 
proportionality  constraints  provides  no  relief  from  the  use  of  target 
values,  then  perhaps  running  the  engagement  as  a  simulation  model,  from 
the  first  day  of  the  engagement,  can  provide  a  more  objective  model  input. 
It  is  this  author's  opinion  that  while  the  use  of  simulation  model  outputs 
for  optimization  inputs  may  be  more  objective,  it  is  doubtful  anyone 
would  vouch  for  the  output  produced.  Further,  the  cost  of  running  such  a 
simulation  in  all  its  required  iterations  might  be  prohibitive. 

Finally,  an  analytical  model  might  be  used  to  provide  input  in  lieu 

of  "subjective"  target  values;  however,  some  significant,  though  more 

subtle,  aspects  of  battle  would  prove  particularly  difficult  to  model. 

In  evaluating  alternatives  for  incorporating  target  activated  munitions 

in  the  NAP  models,  Cudney  and  Bloomquist  state, 

"The  development  of  value  curves  must   be  based  upon  the  judgment 
of  military  commanders  and  analysts  who  are  experiened  in  esti- 
mating the  effects  that  casualties  and  delay   might, have  on  the 
success  or  failure  of  specific  military  missions."   [Ref.  8:  p.  65] 

Combatant  morale  is  also  directly  related  to  choices  regarding  target 
destruction  priorities.  Military  commanders  and  analysts  are  the  sole 
credible  source  for  evaluating  the  effects  of  such  factors. 

Without  consideration  of  such  factors  as  delay  and  combatant  morale, 
this  author  would  find  it  hard  to  accept  any  substitute  for  target 
values  produced  by  an  analytic  or  simulation  model. 
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If  no  apparent,  preferable  alternative  to  the  use  of  target  values 
exists,  why  is  their  use  viewed  so  negatively?  It  is  likely  that  the 
practice  of  manipulating  target  values  to  compensate  for  model  short- 
comings contributes  to  their  poor  reputation.  If  target  values  represent 
just  that,  the  value  of  a  target,  then  military  analysts  would  find  their 
use  easier  to  accept  as  valid.   It  is  difficult  to  stand  behind  model 
output  when  inputs  must  be  manipulated  so  violently. 

Still,  the  issue  of  subjective  model  input  remains.   It  is  this 
author's  belief  the  use  of  such  input  is  allowable,  even  necessitated, 
because  the  use  of  military  judgment  is  still  the  preferred  alternative 
when  trying  to  account  for  all  the  imponderables  existing  on  a  battlefield 
Unfortunately,  while  admitting  military  science  is  a  "soft  science," 
analysts  still  seek  to  "exorcise"  their  models  of  subjective  military 
inputs  (judgments)  rather  than  incorporating  and  exploiting  their  use. 
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IV.   HEAVY  ATTACK  OPTIMIZATION  ALGORITHMS  AND  THEIR  PERFORMANCE 

A.  THE  ORIGINAL  (JUNE  1982)  NONLINEAR  PROGRAMMING  ALGORITHM 

The  optimization  code  present  in  the  June  1982  HA  program  (consisting 

of  those  programs  strictly  subordinate  to  subprogram  OPTMUM  in  Appendix  D) 

appear  to  be  a  faithful  translation  from  the  original  FORTRAN  (to  PL1 

and  back  to  FORTRAN)  of  the  optimizing  code  described  in  [Ref.  2].  That 

code  solves  a  general  nonlinear  program  by  a  sequence  of  local  linear 

programs.  The  algorithm  employed, 

".  .  .  is  a  'local,'  'gradient,'  'stepwise'  correction  descent 
algorithm  ....  By  a  'stepwise'  procedure  we  mean  that  given  a 
point  y  in  the  domain  of  the  functions,  a  'correction'  vector  Ay 
is  determined  and  a  new  point  y  =  y  +  k  ay  is  used  for  the  suc- 
cessor 'step.'   It  is  a  'local'  method  because  the  correction 
direction  Ay  and  its  length  (determined  by  the  skalar  k)  are 
obtained  from  the  behavior  of  the  system  in  a  'sufficiently'  small 
neighborhood  of  the  current  point  y  .  It  is  a  'gradient'  technique 
inasmuch  as  the  gradients  of  the  function  gi(y)  are  principally 
used  to  obtain  the  correction  direction.  .  .  ."  [Ref.  2:  pp.  13-14] 

Use  of  this  algorithm  to  perform  an  HA  optimization  process,  while 
precise  and  correct,  involves  significant  expenditure  of  computing 
resources  (typically  five  to  s'x  CPU  minutes  on  an  IBM  3022).  Anticipat- 
ing a  state-of-the-art  optimizer  would  perform  the  optimization  at  less 
expense,  the  Air  Force  requested  [Ref.  1]  the  June  1982  HA  optimizer  be 
replaced  with  the  X-System. 

B.  THE  X-SYSTEM 

Like  the  original  optimize'",  the  X-System  solves  a  nonlinear  problem 
with  a  sequence  of  local  linear  programs.   (The  HA  internal  model  consists 
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The  objective  function 
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of  a  nonlinear  objective  function  subject  to  entirely  linear  constraints 
so  that  for  each  linear  program,  all  that  is  required  is  that  a  linear 
approximation  be  made  of  the  objective  function  and  then  a  "standard" 
linear  program  is  run.) 

The  X-System  has  been  operational,  but  under  continuing  development 
since  1974  [Ref.  9].  It  is  a  general -purpose,  state-of-the-art  optimiza- 
tion system  which  is  used  as  both  a  vehicle  for  research  and  as  the  basis 
for  a  number  of  commercially  installed,  customized  applications  optimizers. 
It  consists  of  open  FORTRAN  subroutines  and  is  implemented  in  FORTRAN  IV. 
The  subset  of  FORTRAN  with  which  it  is  coded  is  accepted  by  a  majority  of 
FORTRAN  compilers. 

The  X-System  is  designed  to  solve  large-scale  optimization  problems, 
and  is  especially  effective  on  mixed  integer  problems.  Decomposition 
issues  have  been  an  area  of  major  interest  to  the  designers;  however,  the 
core  linear  programming  module  has  received  the  most  design  effort.  It 
exhibits  many  unique  features  including: 

1.  Hyper-sparse  data  representation  [Ref.  10]; 

2.  Complete,  constructive  degeneracy  resolution  [Ref.  11]; 

3.  Basis  factorization  [Ref.  12];  and 

4.  Elastic  range  constraints  [Ref.  9]. 

In  order  to  best  support  the  wide  variety  of  applications  using  the 
X-System,  the  system  is  designed  to  support  all  other  optimization  features 
simultaneously  with  the  nonlinear  feature  (e.g.,  sortie  constraints  form 
an  intrinsic  Generalized  Upper  Bound,  or  GUB,  set  which  has  been  exploited 
in  the  nonlinear  solution). 
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C.   NUMERICAL  EXPERIENCE  WITH  THE  X-SYSTEM  IN  HA 

The  X-System  was  imbedded  in  HA  as  a  subprogram,  without  changing 
the  internal  or  external  models,  and  delivered  to  the  Air  Force  in 
July  1982.  Headquarters,  United  States  Air  Force  (XOX/FM)  has  been 
testing  this  version  of  the  HA  code  (called  "Fast"  Attack  by  its  users) 
since  then.  The  X-System  returns  solutions  whose  objective  function 
values  agree  with  those  provided  by  the  June  1982  code  to  the  second  or 
third  significant  digit  in  less   than  one  fiftieth  of  the  time   while  using 
a  fraction  of  the  compute  region.  However,  the  numbers  of  particular 
sorties  chosen  differ  significantly  in  some  instances.  These  differences 
result  from  the  level  of  precision  specified  for  the  optimal  objective 
function  value  which  permits  early  termination  of  the  algorithm  with  an 
acceptable  solution. 

Since  both  the  optimizers  described  are  supposed  to  deliver  correct 
solutions,  one  can  hardly  help  but  be  surprised  when  their  optimal  solu- 
tions for  the  same  problem  differ,  no  matter  how  small  the  difference. 
One  fundamental  difference  in  the  two  algorithms  is  that  while  the 
original  optimizer  included  a  coded  gradient  function  for  the  objective 
function,  the  X-System  uses  an  automatic  numerical  difference  approxima- 
tion to  estimate  gradients.  Use  of  this  approximation  in  the  X-System  to 
enhance  robustness,  that  is,  to  eliminate  the  errors  and  frustration 
associated  with  coding  derivative  functions,  appears  to  be  responsible  for 
the  small  differences  in  objective  function  value. 

A  solution  precision  factor  is  used  by  the  X-System.  This  algorithm 
parameter  directs  the  X-System  to  stop  optimization  at  the  first  point 
where  the  tolerance  of  the  solution  is  estimated  to  be  comparable  to  user 
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requirements,  or  confidence  in  the  precision  of  input  data.   (This  avoids 
wasting  computing  resources  extracting  the  "last  few  pennies"  in  optimal ity 
from  a  problem  for  which  input  coefficients  have  been  rounded  to  the 
next  higher  dollar.)  By  tightening  this  factor  somewhat,  perceived 
instability  in  total  numbers  of  individual  sorties  can  be  eliminated. 
However,  if  this  tolerance  factor  is  out  of  proportion  with  user  confidence 
in  input  data,  any  "stability"  achieved  is  illusory. 

The  X-System  has  enabled  HA  to  be  used  on  a  time-sharing   system. 
Internal  models  with  81  constraints  (13  GUB)  and  793  variables  typically 
yield  solutions  in  10-15  CPU  seconds  on  an  IBM  3033  processor,  using 
approximately  250K  bytes.  The  much  faster  response  of  the  enhanced 
system  has  not  only  enabled  speedier  evaluations  for  the  theatres,  but 
has  encouraged  a  critical  review  of  the  HA  model.   It  is  now  technically 
possible  to  add  model  enhancements  because  neither  space  nor  time  constrain 
the  problem  as  they  did  with  the  prior  optimizer. 
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V.  A  NEW  HEAVY  ATTACK  EXTERNAL  MODEL  (WEATHER  WARS) 

In  an  attempt  to  ameliorate  some  of  the  major  problems  in  the  HA 
model,  a  new  external  process  was  designed  in  July  and  partially  imple- 
mented in  August  1982.  It  is  still  undergoing  testing  and  refinement. 

The  new  model,  known  as  "Weather  Wars"  (WW),  does  away  with  the  use 
of  weighted  average  sorties  and  attempts  to  logistical ly  constrain  the 
problem. 

WW  avoids  the  use  of  weighted  averages  by  building  sorties  for  a 
particular  weather  type.  Optimization  is  performed  and  inputs  are 
updated  and  reinitialized  as  in  the  former  HA;  however,  at  the  end  of 
processing  for  the  last  time  period,  another  "weather  war"  is  "fought" 
in  which  sorties  are  built  for  a  different  weather  type. 

Running  separate  "wars"  for  each  anticipated  weather  type  (typically 
six)  is  costly.  However,  the  expense  is  offset  by  the  speed  with  which 
the  X-System  provides  internal  model  solutions.  When  weighted  averaqing 
of  weapons  expended  is  performed  following   completion  of  all  "weather 
war"  processing,  the  resulting  weighted  average   solution  is  provided  and, 
at  no  additional  cost,  so  is  the  maximax   solution,  the  number  of  weapons 
required  to  meet  all  constraints  and  inflict  great  damage  upon  the  enemy 
even  if  the  weather  becomes  the  enemy's  consistent  ally. 

WW  logistically  constrains  the  problem  by  noting  maximum  procurable 
quantities  for  each  weapon  type  and  by  reading  in  a  user-defined  number 
of  "best"  weapon  types  for  an  aircraft  operating  against  a  particular 
target  type  in  a  given  weather  condition.  Prior  to  each  optimization, 
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sorties  are  built  and  probabilities  of  effectiveness  (P--'s)  are 
assigned  using  the  "best"  weapon  type  still  available.  If  all  the  "best" 
weapon  types  for  a  particular  aircraft/target  combination  are  no  longer 
available,  the  model  builds  sorties  using  a  fictitious  weapon  type  with  a 
zero  coefficient  of  effectiveness. 

The  internal  model  still  has  no  explicit  logistical  constraint,  and 
given  time  periods  of  fourteen  to  thirty  days,  will  include  sortie  totals 
in  some  time  period  solutions  which  consume  weapons  in  quantities  greater 
than  will  be  available.  However,  this  "over-use"  of  a  weapon  type  can  be 
made  arbitrarily  small  by  reducing  time  period  duration  (and  subsequently 
increasing  total  numbers  of  optimizations).  Another  possible  method  for 
strictly  limiting  the  numbers  of  a  weapon  type  used  is  to  run  the  entire 
WW  model  iteratively,  reducing  the  maximum  procurable  quantity  of  a 
weapon  type  to  the  number  used  in  the  time  period  prior  to  the  one  in 
which  the  "over-use"  occurred. 

WW  permits  analysis  of  realistic  scenarios  in  which  sorties  are  to 
be  selected  after   weather  conditions  are  known.  The  June  1982  HA  infers 
an  assumption  that  sorties  must  be  planned  with  only  synoptic  weather 
forecasts,  or  that  various  weather  conditions  will  exist  throughout  the 
theatre  in  specified  proportions.  WW  presents  opportunities  for  decision 
analysis  incorporating  weather-dependent   target  values  and  damage 
confirmabil ity. 
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VI.   POTENTIAL  FUTURE  DIRECTIONS  FOR  HEAVY  ATTACK 

A.  HA  INTERNAL  MODEL  AS  A  GENERALIZED  NETWORK 

It  is  often  advantageous  to  reformulate  an  optimization  problem  into 
an  equivalent  model  which  may  be  easier  to  solve.  Recent  computational 
advances  in  the  efficient  solution  of  generalized  networks  [Refs.  13,  14] 
have  allowed  these  relatively  specialized  linear  programs  to  be  solved  in 
a  fraction  of  the  time  required  to  solve  them  with  "standard"  linear 
programming  techniques. 

The  internal  HA  model  can  be  viewed  as  a  generalized  network  if 

flight  composition  constraints  are  ignored  and  the  objective  function  is 

simplified  (by  letting  C.  =  0  for  all  j,  or  assuming  a  target  rich 

environment)  as  follows: 

J     I 
maximize:.   £    E   V  . P  .  . S  -  - . 
si(j     J-l   1=1   J  1J  1J 

The  model  can  then  be  viewed  as  displayed  in  Figure  1. 

More  general  views  of  HA,  including  weapons  and  time  periods,  can  be 
formulated  as  multicommodity  compositions  of  generalized  networks. 

The  preceding  perspective  of  the  HA  problem  as  a  generalized  network 
is  of  more  than  purely  academic  interest.  Work  has  been  published 
[Ref.  15]  and  research  continues  [Ref.  16]  solving  generalized  networks 
with  side  constraints  (such  as  flight  composition  constraints).  A 
commercial  quality  optimization  system  exploiting  this  new  technology 
could  provide  a  new  vehicle  for  enhancing  the  efficiency  with  which  HA  is 
processed. 
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Figure  1.  Generalized  Network  View  of  the  Simplified  HA  Internal  Model 


B.   AN  EXPANSION  OF  THE  HA  INTERNAL  MODEL  TO  INCLUDE  WEAPONS  AND  TIME 
Headquarters  Armament  Division  (AFSC/XR)  proposed  to  Headquarters 
United  States  Air  Force  in  January  1982  that  HA  be  reformulated  as  a 
nonlinear  mixed  integer  problem.  Appendix  G  presents  the  basic  concept 
for  the  reformulation.  As  reported  by  Dean  [Ref.  17:  pp.  56,  59],  the 
X-System  has  solved  nonlinear  integer  and  mixed  integer  programs. 
(However,  the  size  of  the  integer  problems  Dean  reported  were  considerably 
smaller  than  the  typical  HA  internal  model.)  The  proposed  formulation  is 
basically  an  expansion  of  the  dimension  of  the  internal  model  to  include 
weapons. 

Unfortunately,  the  model  provided  in  Appendix  G  is  intractable 
as  stated  and  does  not  address  all  the  issues  pertinent  to  HA.   It 
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entirely  overlooks  two  model  components  described  previously:  weather 
and  optimization  over  time.  Also,  it  cannot  ensure  weapon  selection  of  a 
"best"  weapon  for  a  particular  aircraft/target  combination  prior  to  a 
"second  best"  weapon  and  so  on  to  the  "n   best."  A  prioritization 
constraint  is  presented  composed,  in  part,  of  binary  indicator  variables. 
However,  the  mechanism  for  enforcing  the  use  of  these  binary  variables  is 
not  presented. 

While  a  device  for  enforcing  the  use  of  the  binary  variables  can  be 
stated,  this  problem  may  be  solved  at  less  expense  by  allowing  the 
selection  of  sorties  to  occur  over  the  entire  range  of  weapons  considered 
by  the  problem.  If  the  selection  priority  is  ordered  by  P...  (expected 
type  j  target  kills  by  a  type  i  aircraft  loaded  with  type  k  weapon), 
prioritization  is  then  enforced  intrinsically  by  the  presence  of  P... 
in  the  objective  function.  This  alleviates  the  need  for  both  prioritiza- 
tion constraints  and  binary  variables,  so  that  by  expanding  the  size  of 
the  problem,  the  use  of  binary  variables  can  be  eliminated  and  the 
problem  is  made  much  easier  to  solve. 

Using  the  AFXR/SC  proposal,  without  integer  variables,  as  basis  for  a 
new  approach,  optimization  over  time  could  be  addressed  by  further 
expanding  the  dimension  of  the  problem  to  include  time.  The  resulting 
formulation  would  be: 


V-  (K.  ) 
jnv  jny 


J 

i 

ximize: 

£ 

Z 

i  1/  n 

j-i 

n=l 

subject  to: 
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1.     Sorties  available. 


n  NW      /    J  n-1 

n 


<  (1  +  r)   £    RS. 
ft       ln 


i  =  1,    ...,   I;   n  =  1,    ...,  N; 


2.  Target  constraints. 
I    NW 


B-jn  <     .^       SVijkn<Rjn.       J-1.   ...J;n-1.   ....  M; 

3.  Weapon  logistics. 

I     J     N 

£   E   £  Lik  Siikn<Wk>   k  =  l>    •—  NW' 
i=l   j=l   n=l  1K  1JKn    K 

4.  Flight  composition  criteria. 

k  =  1,2,  ...,  k;  m  s  1,  2,  ...  M; 

5.  Target  losses. 

Kjn  =  Dj(n  +  I)' 
All  parameters  retain  their  previous  definitions  unless  (re)defined 

in  the  following: 

T.  -C.        I    NW 

Kjn   =<f  (1  "  *XP  (ljj  UJ"  *  S   &  P1JknV»: 

k  =  weapon  type  index  (k  =  1,  ...,  NW); 

n  =  time  period  index  (n  =  1,  ...,  N); 

r  =  sortie  equality  range  restriction  factor; 

W^  =  number  of  available  weapons  of  type  k; 

RS.   =  resupply  quantity  (initial  quantity,  if  n  =  1)  of  sorties 
available  for  aircraft  i  in  period  n; 
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L-.   =  number  of  weapons  of  type  k  loaded  on  aircraft  of  type  i; 

SijkO  =  °'0;  and 

F...   =  friendly  aircraft  attrition  factor, 
ljkn        J 

Time  periods  would  be  connected  by  a  friendly  aircraft  attrition 
factor,  F...  .  associated  with  each  aircraft/target/weapon  combination. 
The  target  losses  constraints  equate  dead  targets  in  period  n  +  1  (D.,  ,,) 
with  those  killed  by  the  end  of  period  n  (K.  ). 

Unfortunately,  if  any  C.  }   0,  the  target  losses  constraint  is 
nonlinear  and  the  model  becomes  somewhat  more  difficult  to  solve. 
Also,  the  attrition  factor,  F  .    might  realistically  be  defined  as  a 
nonlinear  function  of  time  period  or  prior  attrition.  However,  it  is 
believed  these  nonl inearities  can  be  accommodated  by  a  state-of-the-art 
optimizer. 

The  Weapons  Logistic  constraint  is  appropriate  when  HA  is  used  as  a 
munitions  procurement  model.   In  HA's  secondary  role,  providing  theatre 
commanders  with  insights  into  scenarios,  a  weapon  logistics  resupply 
constraint  might  be  preferred.  Let  W.   represent  the  resupply   of  weapon 
type  k  arriving  in  the  theatre  at  the  beginning  of  time  period  n.  Such  a 
constraint  can  be  stated  as  follows: 
3.1  Weapon  Logistics  resupply. 

I     J 
n  <  CN   *"   i  =  l   j; 


Z        w   -  £   Z         Z    hks.   >0, 

<  CN   <n   i=l   i=l   n  <  CN  1K  JKn 


CN  =  1 N;  k  *  1,  ...,  NW, 

However,  this  constraint,  if  used  as  stated,  produces  a  dense  problem 
matrix.  Model  clarity  and  optimization  performance  would  improve  with 
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the  addition  of  variables  u.      to  represent  the  unused  weapons  of  type  k 

remaining  at  the  end  of  time  period  n.  This  new  constraint  can  be  stated 

as  follows: 

3.1.2  Weapon  logistics  resupply. 

I    J 
\(n-l)  +  Wkn  "  fa     fa   LikSijkn  =  ukn' 

n  =  1,  ... ,  N;  k  =  1,  .. .,  NW; 

where, 

«k0  =  0.0. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  research  was  undertaken  with  the  belief  that  HA  was  a  single 
optimization  model.  Two  models  were  identified  and  characterized.  An 
internal  (optimization)  model  was  found,  and  an  accompanying  external 
model  was  discovered  to  principally  constitute  a  circumvention  of  internal 
model  shortcomings.  Both  models  have  been  analyzed  in  detail,  with  special 
emphasis  on  the  context  of  their  use:  scenario  evaluation  for  theatres. 
The  potential  strengths  and  weaknesses  of  HA  have  been  investigated. 

HA  has  been  provided  with  a  new,  fast  optimizer.  A  new  prototype 
external  model,  Weather  Wars,  which  corrects  some  of  the  more  obvious 
shortcomings  of  HA,  has  been  designed  and  coded.  To  a  degree,  the 
prototype  suffers  from  the  same  lack  of  standard  calling  disciplines, 
variable  naming  conventions,  and  documentation  as  its  predecessor, 
precisely  because  it  is  a  modification  of  its  predecessor.  Any  futher 
attempts  to  modify  the  existing  external  model  are  unlikely  to  provide 
continuing  user  satisfaction.  However,  an  entirely  new  HA  can  be  written 
to  exploit  the  insights  gained  in  this  research. 

What  HA  has  always  needed  is  a  single  model,  a  single  formulation, 
which  encompasses  enough  of  the  real  problem  to  be  adjudged  realistic. 
While  this  may  have  been  impossible  in  1974,  given  the  state-of-the-art 
in  optimization,  it  is  feasible  now.  A  single  model  such  as  that  proposed 
in  the  preceding  section  would  not  only  be  easier  to  understand,  but  when 
implemented  should  actually  run  faster  than  the  current  HA,  providing 
better  solutions. 
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APPENDIX  A.   LINEAR  RESTATEMENT  OF  TARGET  CONSTRAINTS 


Target  constraints  are  of  the  form: 


*j  <  Kj  1  Tj»   J  =  1.  2'  •••>  J; 


which  is  restated: 

J  =  1>  2,  ...,  J. 
First  the  upper  inequality  is  solved  for  the  linear  term  z  P..S.. 
T.  -C.       I 

-C.      I 

(1  -  exp  (  -JL  [a    +  £  P,,S..)))  <  C., 

1  -j   J   i  =  ]_    'J  'J       J 

-C 


ii«"-Tf<vSW»- 


1  -  c 

J 


-c 


iog(i-c.)<(-Ti(aj  +  iEiPijsi.)), 


-T.  I 

J. 
'j 

"TJL  l 

Similarly,  the  lower  inequality  is  solved 
T.  -C.       I 


^nog(i-cj))>(aj  +  ii;iPijsij), 

^OogU-c.))-,.^  P..Sij. 


l.  <  -i(i  -  exp  (  -J  (*  +  £   ?USU)))9 

j   Lj         'j   J   i=i  ^  1J 


c.        -c. 


jTi<(1.exp(^(aj+S  P..S..),), 
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1_iji>(exp(:Ji(.j+  £  P..S..))) 


J 
-c 


J      lj  I  j  J  i=1        1J     U 
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APPENDIX  B.   ALGEBRAIC  MANIPULATION  OF  K.  -  D, 

J        J 

If 

T.  -C.       I 

K.  =  J-   (1  -  exp  (  -J.  (a.  +  £  P  i .-  S  i  . ) ) ) , 

then  noting 

-C.         C. 

exp  (  -J- a  )  =  1  -y^D., 

J   J        J  J 

and  letting 

xo  "  &  pijV 

the  term  K.  -  D.  (from  the  current  HA  internal  model  objective  function) 
can  be  restated  as  follows: 

Kj-Dj=Ji(l-exp(^i(aj  +  Xj)))  -Dj 

T.         C.         -C.  C.    T. 

=  J.  (1  -  (1  -  JL  D.)(exp  (  -J-  X  )))  -  D,  {  J-  )(  J-  ) 

Lj         'j  J        'j  J      J   'j    Lj 

T  C  -C  C 

=  J.  (1  -  (1  -  J.  D,)(exp  (  -J-   X  ))  -  D.(  -1  )) 

Cj         Tj  J        Tj  J     J   'j 

T         C         C  -C 

=  J.  ((1  -  D  JL  )  -  (1  -  J.  D  )(exp  (X.  -J-  )) 
Lj       J  'j         J  J       J   'j 

T.        C.      T.  -C. 

.  J.(l  -  D.  J-)  -  (  J--  D.)(exp  (X.  -J-)) 
Lj      J  'j      Lj    J       J   'j 

■(Jj.OJ).(Jl.DJ)(«xp(XJ^l)) 
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T.  -C. 

=  (  J--  D.)(l  -  exp  (X.  -J-  )). 

Lj    J  J   J 


If  C  has  a  value  of  one,  then  K.  -  D-  simplifies  further  to: 
If  C.  has  a  value  approaching  zero,  then  K-  -  D.  approaches  X-. 

J  vJ  J  J 

Recalling  (1  -  exp  (-x))  is  approximately  x  for  small  x,  and  noting 

-C. 
(  -y.  %.)   ->  0  as  C-   ->  0, 

•J 

and  further  noting  that  as  C-  approaches  zero,  the  quotient  T./C- 
will  become  so  large  as  to  leave  the  effect  of  D.  negligible,  K •  -  D- 

J  J       J 

then  becomes 

VDo  ■ttjH^i»-xj'*«BCj"0- 

In  a  target  rich  environment,  no  matter  what  the  value  of  C,  X-/T. 
will  be  quite  small.  Noting  the  fact  that  (1  -  exp  (-x))  is  approximately 
x  for  small  x,  K.  -  D.  is  finally  approximated  by: 

X.  X.D. 

kj  -  °j  '(Tj-Dj»T7  =  xj-4f 

X. 
=  X.,   for  smal  1  y^-  . 

J  j 

Thus,  K-  is  nearly  linear   in  a  target  rich  environment. 
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APPENDIX  C.     ORIGINAL  HEAVY  ATTACK   INPUT-PROCESS-OUTPUT 


SELECT 
OPTIMUM    .VEAPON 

FOR   EACH 
TARGET.  AIRCRAFT 

COMBINATION 


INDEX   CF 
RELATIVE   TARGET 
IMPORTANCE  IN 

INTERVAL  T 


ALLOCATE   SORTIES   TO 

VARIOUS    TARGETS 

FOR    MAXIMIZING 

MILITARY    .VORTH    OF 

TARGETS   KILLED 


UNIT  PROCUREMENT 

COSTS   FOR   EACH 

.vEAPCN   TVPE 


CALCULATE 

MUNITIONS 

REQUIREMENTS 

FOR    INTERVAL  T 


Inpur/  Ourput 
'actor 


MUNITIONS    MIX 


TARGETS   KILLED 
IN    INTERVAL    T 
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APPENDIX  0.   HIERARCHY  CHART  OF  HA  PROGRAMS  IN  USE  AS  OF  1  JUNE  1982 
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APPENDIX  E.   VARIABLE  NAME  CHANGES  IN  HA  PRIOR  TO  OPTIMIZATION 


Formulat- 

on 

Subroutine 

Subroutine 

Main 

V  ar  i  ab  1  e 

NONLIN 

OPTMUM 

Program 

Si 

SMAX 

S 

SORT 

I 

XNAT 

NAC 

NAC 

J 

XNTT 

NTGTS 

NTGTS 

tj 

T 

TTGT 



°j 

D 

DTGT 

TD 

J 

SL 

TL 



cj 

CC 

CC 

CC 

vo 

VT 

VT 

VT 
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APPENDIX  F.   NAP  MODELS  INPUT-PROCESS-OUTPUT  RELATIONSHIPS 


File  1 
Weapon  delivery  profile  data 


File  2 
Weapon  effectiveness 


File  3 
Payloads 


File  7 
Fractional  kills  per  pass  (FKP) 


File  5 
Target  description 


File  4 
Fractional  kills  per  pass  (FKP) 
by  weapon/delivery  condition 
for  each  aircraft  and  target 
in-year  weapon 


File  3 
Payloads 


y 

FTTe~7 

FKPs  for  each  aircraft  and 
weapon  for  each  target 


File  8 
Attrition  for  each  profile  and 
aircraft  against  every  target 


File  10 
Weapon  profile  and  bomb  damage 
assessment  data 


File  4 
Acquisition  probabilities 


WEAPONICAL 

Computes  the  fractional  kills 
per  pass  (FKP)  for  every 
aircraft  delivery  condition 
and  target  combination. 


INBETWEEN 

Reformates  FKP  data  into  one 
"^  file  with  the  FKPs  grouped  by 
targets 


SURVIVOR 

Computes  the  expected  kills 
per  sortie  (EKS)  and  the 
attrition  sustained  by  attack- 
ing each  target  by  a  specified 
weapon  and  delivery  profile 
with  a  specified  aircraft 
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r 


w& 


e  2 

EKS  and  attrition  for  each 
aircraft  and  weapon  by  target 


File  5 
Cost  data.  Weather  states. 
Delivery  conditions. 


File  8 
Preferred  weapons  by  delivery 
band  for  each  target 


File  5 
Number  of  sorties.  Target 
values.  Constraints. 


SELECTOR 

Selects  the  preferred  delivery 
condition/weapon  for  each  air- 
craft/target/delivery band 
combination. 


Rank  ordering 
of  weapons 


HEAVY  ATTACK 

Identifies  the  sortie  alloca- 
tion (aircraft,  weapons,  and 
targets)  that  maximizes 
mil itary  worth 


|  Sortie,  target  and  aircraft 
I  allocations  that  maximize 
1  mil itary  worth 
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APPENDIX  G.   AFSC/XR  PROPOSED  NEW  HA  METHODOLOGY 
Definition  of  Variables: 

i  =  Aircraft  type  index  (i  =  1,  ...,  I); 

j  =  Target  type  index  (j  =  1,  ...,  J); 

k  =  Weapon  type  index  (k  =  1,  ...,  K); 

T-  =  Number  of  j  type  targets; 

J 

D-      =  Cumulative  number  of  targets  killed  in  prior  time  periods; 

P- -k    =  Expected  number  of  type  j  targets  killed  per  type  i  aircraft 
J      loaded  with  type  k  weapon  when  no  other  targets  of  type  j  have 
been  previously  killed  and  when  conditions  of  conf irmabil ity 
are  perfect; 

S-..    =  Number  of  sorties  of  type  i  aircraft  flown  with  weapon 
J      type  k  against  target  type  j; 

F.     =  Quantity  of  sorties  (fragable)  available  for  aircraft  type  i; 

V.      =  Value  (Military  Worth)  of  target  type  j; 

J 

£ •      =  Lower  bound  on  targets  of  type  j  to  be  killed; 
Q-      =  Number  of  type  j  targets  killed; 

J 

C-      =  Target  kill  conf irmabil ity  parameter  for  j  type  targets 
(as  defined  in  old  formulation); 

W^     =  Quantity  of  type  k  weapons  available; 

L--^    =  Standard  loadout  of  weapon  type  k  on  aircraft  type  i  used 
against  target  type  j; 

n      =  Number  of  members  in  the  weapon  prioritization  sets; 

ri/.  *\      =  The  p   priority  member  (a  k  value)  of  the  ordered  weapon 
p   priority  set  of  k  values  valid  for  the  i   aircraft  type, 
j   target  type  combination  (p  ranges  from  1  to  m   best); 

m      =  flight  composition  constraint  index  (m  =  1,  2,  . . . ,  M) ; 

im     =  aircraft  type; 
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5      =  +1,  maximum  (-1,  minimum); 

J      =  set  of  targets  for  which  a  maximum  (minimum)  flight 
composition  is  required;  and 

9      =  maximum  (minimum)  proportion  of  sorties  flown  by  aircraft 
type  i  against  targets  included  in  set  J 

The  formulation: 

J 


Maximize  Z 

J 


"  S  »J«J  -  DJ 


where 

-C 


Qj  =  cj"    I1  "  ex"  -tJ  <°j  +  S     £pijkS,jk) 


-T.  CD. 


•J  -  Cj  - 

J  ^. 

l"T 

subject  to: 

Sortie  constraint; 

J     K 

E     E  s 

J-l   k-l 

ijk 

<Fi= 

Target  constraints; 

I; 


lj  1  Qj  £  Ty    j  =  l,  ...,  J; 

Side  (flight  composition)  constraints; 


5m 


Ed  -  9  )  s.  .,  +     E(-e j  s1  ik  <  o, 


J~ 


m  =   1,    2,    ...,   M; 
k  =  1,   2,    ...,  K; 


Weapon  constraints; 
I  J 


S      fi  lMSW  £    "k 


;     k  ■  ■  l ,    ...,K; 
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Prioritization  constraints; 


Bi  j  n(i,j)p+1"81  j  n(i,J)p^0' 


1,   •  •  •  »   i» 

j  =  1,  . . . ,  J; 

p  =  1,  .. .,  n; 


where  this  indicator  variable  is  defined  as: 


ljn 


(i,j) 


08  ,f  s«"<u)p "  ° 


p+1    I  1;  otherwise 
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